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ABSTRACT 


The effects of Gaussian-shaped IF and video filters on the 
performance of a square-law detector employing post-detection 
integration are analyzed. The number of additional noise-only 
samples that are integrated due to a finite video bandwidth is 
determined. Emerson’s method is used to obtain an expression 
for the cumulants of the output probability density function. 
These cumulants are used in Edgeworth’s asymptotic series 
expansion of the density functions. By integrating these 
density functions, the receiver operating characteristics are 
determined for various ratios of IF bandwidth to video 
bandwidth. A collapsing ratio which takes into account the 
Gaussian shape of the filters is formulated and compared 
against Barton's approximate formula for the collapsing ratio. 
For typical video bandwidths, Barton's approximate formula is 
found to overestimate the collapsing loss by an amount less 


han 0.5 dB. 
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I. INTRODUCTION 


This report is concerned with the effects of insufficient 
video bandwidth on the detection performance of a square law 
detector with Gaussian-shaped intermediate frequency (IF) and 
video filters in pulsed radar applications. 

The square-law detector is assumed to have an arbitrary 
wide low-pass zonal filter [Ref. 1] so that its instantaneous 
output is just the square of the envelope of the input signal. 
This assumption is in line with Marcum’s definition of the 
detector as "any device whose instantaneous output is a 
function of the envelope of the input wave only" [Ref. 2]. 

In most practical radars with superheterodyne receivers, 
the second detector is usually followed by one or more stages 
of video amplification prior to post-detection integration and 
threshold detection. The bandwidth of the video amplifier 
must be sufficiently wide to avoid distortion of the input 
video signal from the detector. It has been pointed out that 
the video bandwidth B, must be greater than one-half the IF 
bandwidth B,, to be "sufficiently wide" [Ref. 3]. In 
practice, the video bandwidth may be less than half the IF 
bandwidth. 

An insufficient video bandwidth will have two effects on 


the detection statistics of the radar receiver. Firstly, it 


will affect the threshold setting for a specified probability 
of false alarm, P,. Secondly, it will adversely affect the 
input signal-to-noise ratio required to attain a desired 
probability of detection, P,, for a specified probability of 
false alarm. An equivalent statement to the second effect is 
that there will be a degradation in probability of detection 
vis-a-vis the probability of detection attainable with an 
ideal video bandwidth for the same signal-to-noise ratio and 
probability of false alarm. 

The focus of this report will be on the degradation in 
performance associated with the second effect. Mathematically, 
Marcum [Ref. 2] has treated this degradation in performance by 
modelling it as a consequence of integrating M additional 
noise-only samples along with the desired N signal-plus-noise 
Samples’. The degradation in performance is quantified by 
introducing a "collapsing loss", L., which has been defined in 
dB as the incremental increase in signal-to-noise ratio with 
respect to the signal-to-noise ratio required in the ideal 


case of integrating only signal-plus-noise samples for the 


‘Marcum’s paper first introduced the general concept of a 
collapsing loss resulting from the integration of a greater 
number of noise-only samples than signal-plus-noise samples. 
Insufficient video bandwidth is treated, among others, as a 
specific cause of collapsing loss. . 


same probability of detection. Mathematically, this is 


expressed as: [Ref. 4] 


L. dB = Xy,,dB - Xy dB (1) 


where X,,, is the signal-to-noise ratio required when N signal- 
plus-noise samples plus M extra noise samples are integrated, 
and X, is the signal-to-noise ratio required when only N 
signal-plus-noise samples are integrated. 

The signal-plus-noise samples and noise-only samples are 
assumed to be independent and respectively identically 
distributed. 

Mathematically, Marcum showed that for a non-fluctuating 
target, the effect on the detection statistics of integrating 
M additional noise samples can be taken into account by multi- 
plying the joint characteristic function of the N signal-plus- 
noise samples by the joint characteristic function of the M 
noise-only samples. Defining the collapsing ratio, 


M^N 
N (2) 





p =z 


Marcum showed that for the square law detector 





EE 
1+3€ (3) 
(LTS 


exp[ - Np (EM 
Cy, MG) g P 


where C,y..(&) Is the characteristic function for the joint 
density function of N signal-plus-noise samples plus M 
additional noise samples. 

When p is equal to one (i.e., M is equal to zero), Eqn 
(3) reduces to the equation of the characteristic function 
with no collapsing loss, i.e., 


expl -NX( jt ) ] 
a i s 


(EB 





Comparing Eqn (3) with Eqn (4), it is seen that the 


statistical results for p = 1 can be easily extended for any 
p. In particular, receiver operating characteristics (ROC) 
such as Meyer plots“ for p = 1 are available in the 


literature and can be used directly to determine the signal- 
to-noise required by using pN instead of N as the number of 


pulses integrated and with = as a parameter of the ROC 


instead of X.  [Ref. 4] 


^These are plots of P, as a function of the number of 
pulses noncoherently integrated for a given Fehlner’s false 
alarm number and with the input signal-to-noise ratio for a 
single pulse as a running parameter. 


Hence, the quantitative effects on the detection 
statistics of insufficient video bandwidth can be determined 
once the collapsing ratio, p, is known. The problem is in 
relating M (and hence, p) to the video bandwidth. What is 
the effective number of additional noise-only samples that are 
integrated along with the intended number of signal-plus-noise 
Samples given a finite video bandwidth? 

A common practice is to use the following semi-empirical 
formula given by Barton [Ref. 5] to estimate the collapsing 


ratio: 


A 


ae 


PT (5) 


v 


where the IF noise equivalent bandwidth B;; is assumed to be 
matched to the pulse width such that Bt = 1. The video 
filter is modeled as an ideal low-pass filter with a flat 
frequency response from -B, to +B, and zero elsewhere. The 
frequency response characteristics of the actual filter is not 
taken into account apart from computing its noise-equivalent 
bandwidth, B,. 

The purpose of this paper is to investigate the effects of 
Gaussian-shaped IF and video filters on the detection 
performance of the square-law detector. The primary 
consideration in the choice of Gaussian-shaped filters for 


this investigation is the mathematical convenience of the 


Gaussian shape which yields analytical results in closed form. 
Although a Gaussian filter is physically not realizable, 
it has been found to be a good approximation in characterizing 
the frequency response of multistage RC amplifiers [Ref. 6]. 
Hence, the results obtained will have direct application for 
square-law detectors with multistage pre-detection and post- 
detection RC filters, besides providing a relative measure of 
the accuracy of Barton’s formula for the collapsing ratio. 
In the next section, the effects of video bandwidth on the 
detector’s output probability density functions and operating 
characteristics will be analyzed. Specifically, a formula 
will be developed to determine the additional noise-only 
samples integrated due to a finite video bandwidth. This 
leads to the formulation of the collapsing ratio which will be 
discussed and compared against Barton’s formula in Section 
III. Section IV summarizes the methods and findings of this 


paper. 


II. DEGRADATION IN PERFORMANCE 


A. DESCRIPTION OF THE SYSTEM 

The system under consideration is illustrated in Figure 1. 
It consists of a Gaussian IF filter, a square-law detector, a 
video amplifier with a Gaussian frequency response, and a 
linear post-detection integrator. For convenience, it shall 


be referred to as the Gaussian square-law system. 


ri (t) Gaussian t) Square q(t) Gaussian y(t) x y, = z 
IF Law Video E. 
Filter Detector Filter =1 


Figure 1. Block diagram of the Gaussian square-law systen. 


The IF filter has a Gaussian passband characteristic with 
a center frequency of f, and a root-mean-square bandwidth of 
B. The video amplifier also has a Gaussian passband but with 
a root-mean-square bandwidth of v centered at zero frequency. 
Since the gain of the video amplifier will have no impact on 
detection performance, it is assumed to have unity gain for 
mathematical convenience. Correspondingly, the impulse 


response functions of the IF filter and video amplifier are 


respectively defined by: 


1 
Byte) = 2 (208?) * exp[- (2B)? =] cosa,t (6) 


and 


E 
ha (C) EN ZEV E expl- (21v)? 5] (7) 


The signal to be detected is a pulse-modulated sinusoidal 
waveform corresponding to the echo from a target at a 
particular range from the radar. 

The input noise to the system is assumed to be additive 
white Gaussian noise (AWGN) with a two-sided power spectral 
density equal to n/2. The noise-equivalent bandwidth of the 
IF filter, B,,, is matched to the pulse width so that noise at 
the input to the detector can be assumed to be completely 


uncorrelated from pulse to pulse [Ref. 2]. 


B. ADDITIONAL NOISE-ONLY SAMPLES INTEGRATED 

For the hypothetical case of infinite video bandwidth, the 
video range bins will have the same resolution as the IF range 
bins, which is of the order of one pulse width. However, for 
most practical video amplifiers with finite bandwidths, the 
resolution of the video range bins will be poorer than those 
of the IF range bins. In such cases, each sampled output of 


the video amplifier will be shown in the subsequent paragraphs 


to be equivalent to the result of averaging N, independent 
samples of the detector output, where N, is a number greater 
than one. Since the detector follows a square law, each of 
the N, independent samples is equal to the square of the 
voltage envelope in the corresponding IF range bins. 

Suppose a target is present at a particular IF range bin. 
One of the N, samples will consist of signal-plus-noise. The 
remaining N,-1 samples will consist of noise only. Hence, N,-1 
additional noise-only samples will be integrated for every 
signal-plus-noise sample integrated. The added noise results 
in a degradation in performance that is quantified by a 
collapsing loss. 

The number of noise-only samples integrated along with the 
Signal-plus-noise sample can be determined by recognizing that 
the video amplifier with a finite bandwidth behaves as an 
analog filter integrator and paralleling it with a digital 
integrator that averages N, independent samples of the 
detector output. This is illustrated in Figure 2. 

The digital integrator as shown in Figure 2(a) has the 
well-known effect of decreasing the variance of the input 
samples, q,, by a factor of N,, 1.€., 


Var lą,) 
== (8) 


v 


Var ly) = 


where y is the output of the digital integrator. 


Square-Law 


Detector 





(a) Digital Integrator: Var(y) = Var(q)/N, 







Square-Law y(t) 


Video Filter 
Detector 


(b) Video filter: Equivalent Ny = Varíq(t))/Varíy(t)) 


Figure 2. Equivalent number of independent samples averaged 
by the video filter. 


Hence, the equivalent number of samples, N,, averaged by 


the video filter can be determined by evaluating: 


y =. Var ig(t)} 
á Var ly (t)] (9) 
where 
q(t) = output of the square-law detector 


y(t) output of the video filter (amplifier). 


The cumulants of the probability density function (pdf) f, 
at the output of the Gaussian video filter can be determined 
by the method of Emerson  [Ref. 7]. When the input to the 


square law detector is a CW sinusoidal signal contaminated by 


10 


AWGN, the cumulants are given by: 


ONE TAN 


Bu 7 ZZ ———— —— ^ (10) 
MESES (1 + 242, - 1]: 





t 
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1 + za | Lilrzyd)® + aye - DG32y0 ? - 10! (1.5451 
| a EN 
ESL UZ) 1] + 
where 
Ze = ren cumulant of the pdf, f, 
o“ = average noise power at the input to the detector 
X’ = signal-to-noise ratio at the input to the detector 
Bir 
Y 2B 


M 


The mathematical derivation of Eqn (10) is shown in 
Appendix A. For the purpose of this Section, it is sufficient 
to point out that the derivation takes into account the shape 
of the Gaussian filter and that the cumulants are in general 
related to the characteristic function of the associated 


probability density function by the following definition 


meet. 8]: 


exp | b = e | = E (11) 


Zi 


EL 


The cumulants are also related to the central moments of 
the density function.’ In particular, the variance or second 
central moment of the density function is equal to the second 
cumulant of the density function. 


Hence, for the case of a CW sinusoidal signal, 


Var ly(t)) = kR Te mm DR cue | (12) 
vi+2y 


To use Eqn (12) for evaluating N, in Eqn (9), X” must be 
the equivalent CW signal-to-noise ratio that will produce the 
Same variance resulting from the averaging of one signal-plus- 
noise sample and N,-1 noise-only samples. 

Referring to Figure 2(a), let X be the input signal-to- 


noise ratio of the jth range bin containing the target. Then, 


G; (13) 


Mz 


C. I! 


1 1 
= emer m 5 P c 
y Nx Sta 


= 


bu. be 


and 


Var (qj) , i*j (14) 





"A full treatment of the relationship between the 
cumulants and the moments of the probability density function 
can be found in Ref. 8. 
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where g, are independent and identically distributed noise- 
only samples for i = j. 

Now, q, is the output of the square-law detector prior to 
any video filtering. For analysis, this is theoretically 
equal to the output y when the video filter has an infinite 


bandwidth. Hence, 


Varig, - K 


p 2 E EX (15) 
= 02(1+2X) 
and 
Tee we | 
i »,2 | 4-0, x/-0 (16) 


Substituting the above results into Eqn (14) and simplify- 


ing yields the following expression for the variance of y: 


Var {y} = 





ZU NOS 
2 N (17) 
v v 


Let X be the equivalent signal-noise-ratio in each of the 


N, range bins that will produce the same variance resulting 
from averaging the single signal-plus-noise sample with N, 


noise-only samples. 


1 


Then, Using Egn (12) with y =0 And X, it follows 


Var A 
ar tą, Ūsai ) iy 
and 
[d ees 
Var (YU s EIS (19) 
g? = 
= — (1+2X) 


By equating Eqn (17) and Egn (19), it can be easily shown 
that 


N, (20) 


Hence, for an arbitrary y, it follows from Egn (12) and 


Eqn (20) that 


Var WE) DUE. (21) 
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The equivalent number of independent samples N, averaged 


by the video filter can now be determined by using Eqn (9), 





ie.) 
W Var taí(t)! 
V ver) 22) 
with 
ER UE LL. E 
Var ly(t)) We | + - N, 1+y2 S22) 
and 
ee 2x 
var (a(t)? o* (1+ x (24) 


Substituting Eqn (23) and Eqn (24) into Eqn (22) results 


in the following quadratic expression in N,: 


2 
N; + E. t^e. away" din, — 2XV1+2y4 - O (25) 
LY 


Hence, the equivalent number of independent samples 


averaged by the Gaussian-shaped video filter is given by: 


a ee 


; : (26) 


15 


where 


2 
box A ET. (290 
1+y? 
and 
c = -2X/1+2y7 


(28) 


Solving Eqn (26) for y=0 shows that N, is equal to one 
when the video bandwidth is infinitely large. This result is 
consistent with the earlier observation that the video range 
bins will have the same resolution as the IF range bins when 


the video bandwidth is infinite. 


C. RECEIVER OPERATING CHARACTERISTICS 

The collapsing loss caused by insufficient video bandwidth 
can be observed by examining the effect of integrating the 
additional noise-only samples on the output probability 
density function and operating characteristics of the Gaussian 
square law system. 

The number of signal-plus-noise samples N integrated by 
the post-detection integrator corresponds to the number of 
pulses returned from a target as the radar antenna scans 
through its beamwidth  [Ref. 3]. For most practical radars, 
this number is sufficiently large so that the probability 
density function of the voltage at the output of the linear 


post-detection integrator shown in Figure 1 tends towards a 


l6 


normal distribution which can be asymptotically expanded using 


Edgeworth’s groupings of the Gram-Charlier series expansion 


kef. 9, Ref. 10j: 


f (z) = E 





(29) 

Ci (3) / 
= Ž 

a oz) 

C Ga 
T 2 (4) g! + 1 (6) g! 

EX P (2) EE O |z) 

C3 (5) / GC LZ) / ei (9) / 
= => Z — 

wy S DUE 2) serene (2) 
, ] 0 

AS — 

/ÁNK, 2 


where the terms on each line within the brackets are of the 


Same order of magnitude [Ref. 11]; and 


Z = integrator output voltage normalized by the rms 
noise voltage o at the input to the detector. 


lines. Oz - NK, i 





£ 30 
NK, ; 2 
D NEM = (31) 
2m 2 


I7 


Zee 


(32) 
H,(z) = mth order Hermite polynomial in z 
1 K +2 
“= LA E (32) 
N?-5I Goes dies E 


The Symbol RE denotes the rth cumulant of the 
probability density function of the video filter output y 
which takes into account the additional noise-only samples 
averaged along with each signal-plus-noise sample. 

Studies by Cramer [Ref. 11] have validated the asymptotic 
behaviour of Edgeworth’s series and determined that in 
general, the remainder term of the expansion is of the same 
order as the first term neglected. Hence, the output 
probability density function f,(z) can be expanded to the 


desired order of accuracy once the cumulants K are known. 


Yr 


The cumulants K,, can be determined from Eqn (10) with X’ 


Y,r 


X 
equal to rb l.e., 


v 


(20?)* (r-1)! 


ee HS (33) 
[yY1+2y2 + 1]7 - [y1+2y? - 1] * 






EN] 
į (1+ m 
+ 1]? + L(1+2y2) 2 - 1] 7 
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Once the probability density function f,(z) is determined, 


the cumulative distribution function F,(z,) can be found by 


integrating f,(z,) in Eqn (30) from - to z,. Hence, 





(34) 
E S gia (2!) 
T = $9 (z/) + = $9) (z/) 
- QG. ec (zh - ER ato (uh - 900 (2) 
pre] 
where 

E Sum NA nl 
NK, 2 

P(z/) SHE ew - dide 0.5 ert} =| (36) 


The probability of false alarm P,, and the probability of 
detection P, for a normalized threshold voltage of zZ, can then 


be determined by taking the complement of Eqn (34), i.e., 


Pe, E Lo) (37) 


O 
a. 
| 


ee (38) 
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where F,.(2,) denotes the cumulative distribution for the 
hypothesis that the target is absent while F,,(2,) denotes the 
cumulative distribution for the hypothesis that the target is 
present. 

The receiver operating characteristics is then obtained by 
plotting the probability of detection versus the probability 
of false alarm with X, y, and N as parameters. 

It is obvious from the above expressions that the process 
involved in the computation of the output probability density 
function f(z) and the receiver operating characteristics is 
numerically tedious. However, the steps involved are straight 
forward and can be easily programmed for execution in a 
digital computer. As part of this thesis, routines in PC- 
MATLAB have been written to compute the cumulants, the 
probability density functions, and the associated receiver 
operating characteristics. These routines are listed in 
Appendix B. 

Using these routines, the effect of insufficient video 
bandwdith on the the detection statistics of the Gaussian 
square-law system can be investigated for various values of Y, 
input signal-to-noise ratio, and N. This is illustrated in 
Figure 3 to Figure 5. 

Figure 3 compares the output probability density functions 
obtained for a finite video bandwidth equal to a quarter of 


the IF bandwidth (y=2) against the corresponding density 
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Figure 3. Gaussian square law system output probability 
density function. 

functions when the video bandwidth is infinite (y=0). The 
output probability density functions corresponding to the 
absence and presence of a target are respectively denoted by 
f,..(z) and f,,(z2). For this illustration, the number of video 
Samples integrated, N, is chosen to be nine and the signal-to- 
noise ratio in the IF range bin is assumed to be 3 dB when 
there is a target return. The choice of these parameters is 
in part to facilitate comparison with Marcum’s plots [Fig.4, 
Ref. 2] which assumed infinite video bandwidth. Accordingly, 
the plots in Figure 3 for Y-0 are found to be in agreement 


with Marcum’s results which serves as a validation of the 
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method and routines developed to compute the probability 
density functions. Graphically, the effect of the finite 
video bandwidth is observed to be two-fold. Firstly, the 
variance of f,,(Z) is reduced vis-a-vis the density function 
for infinite video bandwidth. Secondly, the variance of 
f,.(2) is also reduced but this is accompanied by a leftward 
shift of the mode. The net result is that for the same 
probability of false alarm, the detection performance has 
degraded considerably. This is the consequence of the 
integration of additional noise-only samples induced by the 
finite bandwidth of the video amplifier. 

Figure 4 shows the receiver operating characteristics for 
y=0 and y=1. The number of video samples integrated N is 
equal to 64. For a given signal-to-noise ratio and 
probability of false alarm, it can be seen that the 
probability of detection achievable for the system with the 
finite video bandwidth (y=1) is poorer than the system with 
the theoretical reference of infinite video bandwidth. 
Conversely, the system with the finite video bandwidth will 
require a higher signal-to-noise ratio to have the same 
performance as the system with infinite video bandwidth. This 
is illustrated in Figure 5. For a probability of false alarm 
and probability of detection equal to 10” and 0.95 
respectively, a system with infinite video bandwidth will 


require an input signal-to-noise ratio of O dB. To achieve 
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Figure 4. Gaussian square-law system operating 
characteristics 
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Figure 5. Gaussian square-law system operating 
characteristics. 
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the same performance, the system with a video bandwidth that 
is half the IF bandwidth will require an input signal-to-noise 
ratio of 0.61 qB. From Eqn (1), this implies a collapsing 
loss of 0.61 dB that is attributable to the ‘insufficient’ 


bandwidth of the video amplifier. 
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III. COLLAPSING RATIO 


The statistical analysis in the preceding section provides 
the theoretical background and insight concerning bandwidth 
effects of a Gaussian-shaped video filter on the performance 
of the square-law detector. The analysis has also shown how 
sets of receiver operating characteristics such as those shown 
in Figure 4 can be generated using a digital computer to 
determine the input signal-to-noise required to attain a 
desired probability of detection for a specified probability 
of false alarm and a given IF bandwidth to video bandwidth 
ratio. While such an approach is feasible, it is more 
convenient in practice to follow Marcum’s approach, using the 
collapsing ratio p defined in Eqn (2) to determine the 
required input signal-to-noise ratio. This is because the use 
of the collapsing ratio will obviate the need to generate sets 
of receiver operating characteristics for each value of y. 
As explained in Section I, only the receiver operating 
characteristics for an infinite video bandwidth (y=0) are 
necessary once the collapsing ratio is known. Furthermore, 
there are usually other sources of collapsing loss apart from 
insufficient video bandwidth in practical systems. In such 
cases, the formulation of the collapsing ratio induced by each 


source will facilitate computation of the overall collapsing 
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loss. A general treatment of the problem of multiple sources 
of collapsing loss can be found in Meyer [Ref. 4] and Barton 
[Ref 5]. It is however beyond the scope of this report. In 
this section, the aim is to formulate the collapsing ratio 
associated with insufficient video bandwidth for the Gaussian 
square-law system and to compare the result with Barton’s 
formula for the collapsing ratio in Eqn (5). 


From Marcum’s definition of the collapsing ratio, 





P > cu (39) 


where M is the number of additional noise-only samples 
integrated along with the N signal-plus-noise samples. 

From the analysis of Section II, the number of additional 
noise-only samples integrated for every signal-plus-noise 
sample integrated in the Gaussian square law system is equal 
to N,-1 where N, is given by Eqn (26). It follows that if the 
number of signal-plus-noise samples integrated is equal to N, 


then 
M = N Se 1.) (40) 
and 


Po = La a N, (41) 
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Hence, the collapsing ratio caused by insufficient video 
bandwidth in the Gaussian square-law system is equal to the 
equivalent number of independent IF samples N, averaged by the 
video filter. 

The subscript ‘G’ has been used to denote that the 


collapsing ratio p, given in Eqn (41) is caused by insuffi- 
cient video bandwidth in a Gaussian square-law system. It is 
also used to distinguish from Barton’s formula for the 
collapsing ratio: 


Bir m 2B, 


ae og SY (42) 


v 


For convenience, p, and p, shall respectively be referred 


as the Gaussian system formula and Barton'/s formula in the 
context of collapsing ratios caused by video bandwidth. As 
pointed out in Section I, Barton’s formula models the actual 
filters with ideal rectangular shaped filters having band- 
widths equal to the noise-equivalent bandwidths of the actual 
filters. A comparison of the two formulae will provide a 
relative measure of the accuracy of Barton’s formula in 
modelling systems that approximate the Gaussian square-law 
system. 


Figure 6 plots the collapsing ratio formulae p, and p, as 


a function of the ratio y of IF bandwidth to twice the video 


bandwidth.  |The plot obtained by applying Barton's formula is 
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Figure 6. Collapsing ratio due to insufficient video 
bandwidth. 


indicated in the figure. As expected from Eqn (42), it is a 
linear plot and is independent of input signal-to-noise ratio. 
All the other curves are obtained by applying the Gaussian 
system formula for input signal-to-noise ratios of 0, 0.1, 1, 


and 2. It is is seen that p, is a function of the input 


signal-to-noise ratio X. For the same y, the collapsing 


ratio given by p, decreases as the signal-to-noise ratio 


increases. This dependency complicates the method of using 
the collapsing ratio to determine the required signal-to-noise 


ratio to attain a specified probability of detection and 
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probability of false alarm since the collapsing ratio is 
itself a function of the signal-to-noise ratio. For very 


small signal-to-noise ratios (X<1), the collapsing ratio p, 


can be approximated as 
Qu = /1 at DE (43) 


which is independent of the input signal-to-noise ratio. 

The collapsing ratio given by Eqn (43) can be considered 
as the worst case situation for any given y in the Gaussian 
square-law system. Taking p.-vi+2y? as the reference, it is 
seen that Barton’s formula overestimates the collapsing ratio 
for y less than two. In this region, the largest relative 
deviation is about 22.5% in the vicinity of y=0.5. For y>2, 


Barton’s formula underestimates the collapsing ratio. For 
large Y, Ps/P¿ approaches 1//2 so that Barton’s formula 


underestimates the collapsing ratio by a factor of about /2 


for very large y. 

The collapsing loss associated with each collapsing ratio 
can be determined by finding the incremental increase in 
signal-to-noise ratio required with respect to the signal-to- 
noise ratio required in the ideal case of infinite video 
bandwidth for the same detection performance. Marcum has 
found that the collapsing loss is a function of the number of 


signal-plus-noise samples integrated but varies only slightly 
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with the probability of detection P, and probability of tatce 
alarm P,  [Ref. 2]. His results for the collapsing loss as 


a function of the collapsing ratio y and with Nas a parameter 


is shown in Figure 7 for P,-0.5 and P, -10710, 


Collapsing loss Le (db) 





1 2 4 40 lOO 300 000 | 
Collapsing ratio p 


Figure 7. Collapsing loss vs. collapsing ratio for 
fons and P. OS) (From Ref. 4) 


Using Figure 7 and the respective equations for p, and p,, 


the collapsing losses predicted by the application of Barton’s 
formula are compared against those deduced using the Gaussian 
system formula for selected values of y in Table I and Table 


II. The collapsing loss associated with the use of Barton’s 
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formula is denoted in the table by L,, while that associated 
with the Gaussian system formula is denoted by L,,. Using L. 
as the reference, the difference between the two collapsing 


loss for a given y is denoted by aL.,. In Table I, the input 


signal-to-noise ratio X is assumed to be much less than one so 


that the collapsing ratio p, can be taken to be equal to 


y1+2y*. A signal-to-noise ratio equal to one is assumed in 


Table II to provide a relative bound on the collapsing loss 
for signal-to-noise ratios between zero and one. 

The comparison shows that for small signal-to-noise 
ratios, the collapsing loss predicted by Barton’s collapsing 
ratio is within 0.3 dB of that deduced by the Gaussian system 
formula. If the signal-to-noise ratio is near or equal to 
one, Barton’s formula overestimates the collapsing loss by 
about 0.4 dB for video bandwidths around half the IF 
bandwidth. 

Typically, the number of pulses integrated in a search 
radar may range from thirty to forty. In such cases, the 
required single-pulse input signal-to-noise ratio is close to 
unity for reliable detection performance [Ref - „air From 
Table II, it follows that the collapsing loss will be 
overestimated by about 0.4 dB if Barton’s formula is used in 
a Gaussian square-law system with a video bandwidth about half 


the IF bandwidth. 
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TABLE I. COMPARISON OF COLLAPSING LOSS. 
BARTON’S FORMULA VS GAUSSIAN SYSTEM FORMULA. 


Collapsing loss, Lc dB 


Input SNR 
X «1 





TABLE II. COMPARISON OF COLLAPSING LOSS. 
BARTON’S FORMULA VS GAUSSIAN SYSTEM FORMULA. 


Input SNR Collapsing loss, Le dB 





IV. CONCLUSIONS 


The effects of Gaussian-shaped IF and video filters on the 
performance of a square-law detector in pulsed radar 
applications employing post-detection integration have been 
analyzed in this paper. 

Specifically, a formula for determining the number of 
additional noise-only samples integrated due to the finite 
bandwidth of the video filter has been derived.  Emerson's 
formulation for the cumulants of the voltage probability 
density function at the output of the video filter has been 
used in this derivation. 

The resultant cumulants of the probability density 
function at the output of the post-detection integrator are 
obtained. Based on these cumulants and Edgeworth’s asymptotic 
expansion of the Gram-Charlier series to represent the density 
functions, MATLAB routines have been written to plot the 
receiver operating characteristics for various ratios of IF 
bandwidth to video bandwidth. The degradation in performance 
or collapsing loss can be observed and quantified from these 
operating characteristics as the video bandwidth is reduced 
relative to the IF video bandwidth. 

Taking into account the Gaussian shape of the filters, 


the collapsing ratio is found to be equal to the equivalent 
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number N, of independent IF samples averaged by the video 
filter. This is given by Eqn (25) which is a quadratic 
equation in N, with signal-to-noise ratio and IF bandwidth to 
video bandwidth ratio as parameters. The collapsing ratios 
given by Eqn (25) and the associated collapsing loss are 
compared against those predicted by Barton’s formula in Table 
I and II. The comparison shows that for typical signal-to- 
noise ratios close to O dB and video bandwidths about half the 
IF bandwidth, Barton’s formula for the collapsing ratio 
results in an overestimation of the collapsing loss by about 


0.4 dB. 
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APPENDIX A 


MATHEMATICAL DERIVATION OF THE CUMULANTS 
This appendix follows the method of Emerson [Ref. 7] in 
deriving a closed form expression for the cumulants of the 
probability density function f,(y) at the output of the 


Gaussian video filter shown in Figure A-1. 
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Video 
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Square 
Law 
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Figure A-1. Block diagram of the receiver. 
The input to the system is a Gaussian random process given 
by: 
mee) E Ss (E) + n(t) (AFL) 
where 


s(t) - /2A coso,t (A-2) 


and n(t) is additive white Gaussian noise (AWGN) with a power 


spectral density equal to m/2. 
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The impulse response functions of the Gaussian-shaped IF 


filter and video filters are respectively defined by: 


i 
h;,(t) mes 2 (2xB-) ^ exem (228)? —-] cos2T fab (A-3) 


and 


i 
12 (27v?) ? exp [- (21v)? ] | (A-4) 


where ß and v are respectively the root-mean-square (rms) 
bandwidths of the IF and video filters. 

The starting point in this derivation is the following 
general expression of the cumulants K, obtained by Emerson for 
an arbitrary input signal contaminated by AWGN into a square- 
law detector with IF and video filters of arbitrary frequency 


response characteristics: 


K, = (qn) AES Sai (n)t tr! So Ags, (t)? (A-5) 
ze i-0 


where 


oo 


s,(t) - [s(t-x) b; (x) dx (A-6) 


= 09 
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Tea e are respectively the eigenvalues and eigen- 


functions defined by the integral equation: 


A;o;(x) = [ $66» $, (y) dy (A-7) 


The function g(u,v) is called the system kernel and it is 


completely defined by the impulse response functions of the 


filters. In particular, it is defined by the following 
equation: 
g(u,v) - f fic (u- 2 A E (A-8) 


- 00 


which can be expanded into the uniformly convergent bilinear 


Mercer series, i.e., 
g(u,v) = > Azo, (u) o; (Vv) (A-9) 
i-0 


To determine the cumulants for the Gaussian system shown 
in Figure A-1, the method of solution is to evaluate the 
system kernel using Eqn (A-8) and expressing it in a form that 
matches the Mercer series. The eigenvalues and eigenfunctions 
that satisfy the integral equation of Eqn (A-7) can then be 
identified. The cumulants can then be determined by using 


these eigenvalues and eigenfunctions to evaluate Eqn (A-5). 
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Thus using Eqn (A-8), the system kernel for the Gaussian 


system shown in Figure A-1 is given by 


g(u,v) = [ 2vax Bexp[- (2nB)*42)" ] cos o, (u-z) 


y Vm vexpl- (22v)? | 


x 2y2mx Dexp[- (22)? EL | cos w, (V-Z) 


pm 


x dz (A-10) 


Now, 


cosw,(u-z) coso,(v-z) 


(A-11) 
=. (u-v) 
= = le tme + 
ee ( Em + Si (u+v) sin2 | 
> sino,(u+v) sin20,Zz + sin, (u+v) sin20,Z | 
Substituting Eqn (A-11) into Eqn (A-10), 
> OU, vie An D? /2n ien 4 e (à-12) 
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where 


J, = COs, u-v) f exp | - mh! Waz) (v-z)?] 
+ Sen 2T dz (A-13) 
go = f coso,(uev-22) exp | - (GP p (u- 2)? (v-2)*] 
Eoi 3 ES (A-14) 


The integral in Eqn (A-13) can be evaluated by completing 


the square in the exponent and casting the integrand into the 


form of the standard integral identity: [Ref. 12] 
f expl-a? (x-m3) dx = vu 
a (A-15) 


-m ` 


Thus completing the square, 


C uva) run) 





= cosw_ (u-v) e - 
Jı HE L 2 1+21 1+2y* 
p 242 2 
x f exp 1-22 É (14239) [z- Nese uv) io) ez 
y? 1+2y? 


— 00 


(A-16) 
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where 





(A-17) 


The symbols B, and B, in the above equation are 
respectively the noise equivalent bandwidths of the Gaussian 


IF and video filters. 


i : 2 D 
Making use of Eqn (A-15) with a = =~ (1+2y°), 
Y 


|  (2nf)? , 1-«y? 2y?uv 
= coso.(u-v) exp | LE a4 ty) - eae 
gi Oo P | 4 1+2y? 1+2y? 4 
x L 
a (A-18) 
/2n v (142wy?) ? 
NOW, 


consider the second term g, 


of the system kernel 
given in Eqn (A-12). 


By completing the square in the exponent 


as in the evaluation of g,, it can be similarly shown that 


2 2 2 
J, = exp poate (rya 
4 1+2y? 


1+2y* j (A-19) 


x 
+ —8 


cosw,(u+v-2z) exp[-a*(z-m)*] dz 








where 
RnB 2 
a = = (Iray JS (A-20) 
Y 
m = CEV 2 
ES (A-2309 


40 


Now, 


f C O a E -a-(zZ-m)-] az (A-22) 
= f cosw,(u+v) cos2w,Z exp [-a?(z-m)?] dz 
+ f sinw,(u+v) sin20,z exp[-a*(z-m)?] dz 
mec x= (z-m). Then, 
(A-23) 


i cosw,(u+v) cos2w,z exp[-a*(z-m)*] dz 


eo 


cos o, (u- v) cos2wm | cos20,x exp[-a*x*] dx 


— 00 


oo 


- cosQ,(u-v) sin2egn | sin20,x exp[-a?x*] dx 


From Eqn (679) of Ref. 12, 





a 2 
| cos20,x exp[-a?^x*^] dx -= YE exp | - = | (A-24) 
= a“ 

(A=25) 


f sin20,x exp [-a*x*] dx = 


- 0 
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Hence, 


| coso,(u«v) cos20,z exp[-a^(z-m)*] dz 





(A-26) 
VE expl... 
= cosw,(utv) cos20,m YE exp | - = 
a a? 
The second integral in Eqn (A-22) 


can be similarly 
evaluated to give the following result: 


o 


[ sine, (Cus) sin20,z exp[-a*(z-m)*] dz 


-09 


(A-27) 
= Sino (u*«v) sin2 yu We Qo 

E 0 „M Exe mi 
Using the results of Eqn (A-22), Eqn (A-26), and 

Eqn (A-27) to evaluate Eqn (A-19), it can be shown that 
| (omn)? 1+y? 2y2uv y | 
> = ex -JANRI E E £ 
2 g 4 1+2y* ED ye | (A 2 
1+y* 
e ee) 
x exp | = 


— À — ee 
2n2ß? (1+2Y2) a 
D "i 2T v (1+2?) ? 


Comparing Eqn (A-28) with Eqn (A-18), it is obvious that 


if w, > 2nB, then g, will be small and negligible compared to 
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Oe. Hence, neglecting the contribution from g,, the system 
kernel given in Eqn (A-12) becomes 
4nB*cosw,(u-v) 


Gon NN e EE (A29) 


(ay 


2 2 2 | 
x ep [ceo (u24y2) ..2Y uv) 
2 1+2y? 1+2y? ? 
Using the following substitutions, 
i 
INPS y A (ESO) 
2n 
= 
ee ey (A-31) 
2np 
1+2y*-1 


(A-32) 
y1+2y? +1 


Eqn (29) can be transformed into the following form suitable 


for the application of Mehler’s formula: 


g(u,v) - [cose,ucoso,v - sino,u sino,v ] (A-33) 
2 2 Ž [ Ž 2 ZN 1 
- An exp | - Ay leol tuy» D 


/142y? 2 1-t£? J 


Using Mehler’s formula, 





2 2 E 3 Ž 
exp nl ae | = y1-t?* 2; H; ( X) A; (y) d 


(A-34) 
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where H,(x) is the Hermite polynomial defined by 


i d? 
From Eqn (A-30) and Eqn (A-31), 
anpu 
= u B (A-36) 
(name 
2nBv 
d Nm (A-37) 
(1:522): 
Substituting the above equations into Eqn (A-33), 
gius v) = Goa Ga (A=38) 
where 
DA 
g,(u,v) ¿ME ze (A-39) 
1+y1+2y? V1+2y2+1 


ls 


ql 
[022 11 (Ta 


| (aro) LZ a 
/AnD exp, --———t——— — 
a 2 | 
[ 224! Ym (1+2y?) 4 ] 
> y H, | 2npu EA 2n pv 


La T H 
(Ira k C ee 


44 


es ee ap (A-40) 
1+y1+2y* VD 


En, S 


| __(2peß)? u? | 
varß exp | A | cose,u 
i 1 
CA a 2. ] + 
EE a L 
i VAnp EXP | ATA > | S1nOo,V 


icu 
A T (t) 


= | 2rßu 1, | 217 Bv 
le [er] 
i (1+2y2) 4 (1*2y?) 4 


By comparing Eqn (A-38) against the Mercer series 


expansion of the system kernel 


ao 


g.(u,v) = $, X,i;6,;Q00 6,100 (A-41) 


i-0 


the eigenvalues and eigenfunctions can be deduced to be 


ME by ay td (A-42) 


y 12 A23) = fl say? +1 


yarß exp | _ mp)? e | H; | nn. | COS QU 
2 = 
$ ME ed (1+2y?) : 


1 


] 2 


E 
4 


Ec TES 


(A-43) 
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Similarly, it can be shown that 


s. i A Àc i (A=) 
f o (22rb) u loo 2nBu |_. 
/AnB exp, -————— — | 4; sinw,u 
A eee 
6. „(u) = y) 
2 i ı 
(211121 
(A-45) 
Hence, from Eqn (A-42) and Eqn (A-44), 
(Ae)? y (Ag (A-46) 
170 1-0 
BYE y p yYırzy? 1 qu 
(1+y1+2y?)* i-o y1+2y* +1 
y 4 
Now, t P gm for MEAN 
Hence, 
= E - 2 (2B0/n) * 
(A. ;)7 = AA ee, og 
2 l 2 (Joye 41 ety 14 oye Sa 
(A-47) 
AA ^" >» Cen CE E 
i-0 i-0 i-0 
pa 2(2B/T) (Aaa 


(/142y? «1) 7- (J1*2y? -1)* 
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From Eqn (A-12) and Eqn (A-45), 


oo 


Se i(t) = /2A [ cosw,(t-u) $. ,(u) du 


—00 


al 


„ZA (4nB) ž 
[es 
[2izl!ym (142y?) * ] ? 


(A-49) 


oo 


2 
x [ exp (-a? 4) H,(a@u) cosw,u cosw,(t-u) du 


— 00 


where 
E 2n 
T (250) 
(1+2y4) 4 


cosw,u coso,(t-u) = COS Wat + = cos (2w, U+ t) 


(A-51) 


For 0o,»2np, the residue from the high frequency 


component in Eqn (A-51) can be neglected in the integration of 
Eqn (A-49). Hence, 


1 
/2A(Anf) ? cos w,t (A-52) 
1 2 


8. ¿ (€) 1 
( 22a: fm (1+2y2) 4 ] 2 


p 2 
X f expl -a= ) H;(«u) du 
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The following integral identities derived in 


useful in evaluating (A-52): 


exp (-x*) Han. (Y2x) dx = 0 


g 8 


exp (-x*) H,,(/2x) dx = Lam lyr 


D ME 


2 
Let aE 


Then au = /2x and du = dx. 


Ref. 


13 are 


( A-9589 


(A-54) 


(A-55) 


Using the above substitution, Eqn (A-53) and Eqn (A-54) 


oo k e 
| exp(-a^—-) H;(«u) du 


- 0 


Io 
XL [ exo (2x) E, (220. dx 


Hp 


LE (1+2y2) % i! 
2rB (2) for i even 


I 


¿ 
| 
| 0 for i oda 


48 


(A-56) 


Hence, using the result of Eqn (A-56) to evaluate Eqn (A-52) 





al 
> E S ; 
2A )? (1«2y2) 8 (1!) cos o,t for i even. 
2B/T a 
27 (=)! 
2 
| 
EN (t) = 
y | (A-57) 
| 0 for i odd. 


A similar result is obtained for s,,(t) except that the 


cosw,t term in Eqn (A-57) is replaced by a sino,t term. 


Using Eqn (A-40) and Eqn (A-55), it can be shown that 


= r- i 
DE (4) sealt) = ar Mau “E20 ) ees. GW) (e 
-œ (1-/1*2y?)7 


go yop A jar GDI 
= o 22i (i1)? 


(A-58) 
Let [ X1*2Y'-1 Jar - X (A-59) 
y1+2y? +1 
From Ref. 12, 
T, 1 
oo f ae a? 
ea) il 2 
E E = — 
i-0 5 222 i 2^7 2! : 
-1 
eme (A-60) 
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Hence, using this result to evaluate Eqn (A-58) gives 


T 
= 2 r-1 2y 4 
y VA 0 see) = A 2 AA EA costw,t 
[-0 


i (WTA 


(A-61) 
Similarly, it can be shown that 
E i 
2 ED 2\ 4 
y (A cT Se ae) = = see) ea sintw,t 
db [ G13237 «1)27- (13237 -1) 27 ] 2 
(A-62) 
Hence, 


eo 


(Az) * si (t) > Mas) ELO ee a CE ES 
i-0 


1-0 


M: 


H 
a 
o 


i 
_ 22° (2ßvVr) il (1+2y7) 4 


a 


[ (y1+2y? +1)27 e (MED? — 3125 ] 2 


(A-63) 
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From Eqn (A-5) and the results of Eqn (A-48) and Eqn (A-63), 


p DAY (mur Y Ap site) 


i-0 i-0 


K, = (mn 


(2nBvT) 7 (r-1)! 
(V1+2y2+1): Tay -1)* 





2 


npyr 


(V1*2y7«1)7- (/1:2y? -1)* 
1 
[(Vi+2yŽ+1)“- (Vl+2y“-1)5 J ? 


Paez 





i 
x (1+2y2) | (A-64) 


Now, it can be shown that the IF noise bandwidth B,, is 
equal to Byr. It follows that the total average power into 


the detector is given by 


c? = nfyx (A-65) 

Hence, 
K = BER car |) eee x A-66 
aay eye = (vy ieoy7=1) 7 ( ) 


| 14 120 | (1:2y *1)' - (V1:2y -1)7 (142?) * | 
o? N (Yırzy?+1)? + (/1+2y2-1)* 
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This equation which expresses the cumulants K, in terms 


of the IF bandwidth to video bandwidth ratio y and the input 
signal-to-noise ratio A*/o* concludes the mathematical 


derivation of this appendix. 


52 


APPENDIX B 
MATLAB PROGRAMS FOR PLOTTING 
THE RECEIVER OPERATING CHARACTERISTICS 

This appendix is a listing of the Matlab programs written 
as part of this thesis to plot the probability density 
functions and operating characteristics of the Gaussian square 
law system described in Section II. The programs are written 
as MATLAB M-Files. The following are respectively the 
filenames of the programs for plotting the probability density 
functions and operating characteristics: 

1. SQ2PDF.M: This plots the probability density functions 
based on Eqn (29) which is an Edgeworth’s asymptotic 
expansion of the Gram-Charlier series. 

2. SQ2ROC.M: This plots the receiver operating characteris- 
tics based on Eqn (37) which is derived by integrating 
Eqn (29) from the threshold voltage to infinity. 

These are script files which are run by invoking their 
filenames in the MATLAB operating environment. Both programs 
will prompt the user to enter the necessary input parameters 
prior to execution. In addition five function type M-Files 


are specifically created to facilitate programming and 
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execution of these script files. The filenames of these five 


functions and their purpose are: 


1. GCUMU.M: This function computes the cumulants of the 
video filter output probability density funtion based on 
Eqn (10). 


2. GDFZ.M: This function computes the probability density 
for given normalized values of the Gaussian square-law 
system output based on Eqn (29). 


3. UPFZ.M: This function computes the probability of false 
alarm or probability of detection of the Gaussian square- 
law system for a specified threshold voltage. The 
computation is based on Eqn (37) and Eqn (38). 


4. GCR.M: This function computes the collapsing ratio 
caused by the finite video bandwidth. It is based on Eqn 
(26). 


5. HMT.M: This function computes the mth order Hermite 
polynomial for values of m less than or equal to 12. 


The MATLAB source programs of the two script files and the 
five function files are provided in the remaining pages of 


this appendix. 


54 


1.  SQ2PDF.M PROGRAM 


This program computes and plots the pdf of a square law 
detector with Gaussian IF and video filters. It will 
prompt the user to enter the following parameters: 

a. Range and resolution of the x-axis for the plot. 

b. SNR = input signal-to-noise ratio. (Not in dB) 

C. GAMMA - IF bandwidth to twice video bandwidth ratio. 
d. N = Number of pulses to be integrated. 

In addition to the pdf plot for the specified GAMMA, the 
pdf plot for an infinite video bandwidth will automatically 
be plotted for comparison purposes. 


oo o9 oo oo o? o9 o? o? o oe 


oe 


disp('The output pdf will be plotted as a function of the’ 
disp(’normalized voltage, 2.') 
disp('Enter range and resolution of z to be plotted’ ) 


range = input(’{zmin,zmax,zres] : ’) 

N = input(’Number of pulses to be integrated =  ') 
SIGMA = input(Input noise rms voltage = D 

SNR = [0 0] 

GAMMA = [0 0] 

ORDER = 4 


SNR(2) = input('Enter input SNR:  ') 
GAMMA(2) = input(’Enter value of GAMMA: ^") 


drsp(l “Computing ...’) 
KYO = [] 
KY1 = [] 


Z = range(1):range(3):range(2) 
GAMMAO = GAMMA(1) 
for j = 1:2 
iK — 
SNRj = SNR(j) 
for r = 1:ORDER 
K(r) = gcumu(r,SIGMA,SNRj,GAMMAO) 


end 
if j == 1 
fZinfO = gdpz(ORDER,N,Z,SIGMA,SNRj,GAMMAO) 
else 
fZinf1 = gdpz(ORDER,N,Z,SIGMA,SNRj,GAMMAO) 
end 
end 
GAMMAv = GAMMA(2) 
SNR} = 0.0 


for r = 1:ORDER+2 
KYO(r) = gcumu(r,SIGMA,SNR7,GAMMAv) 
end 
fZvo = gdpz(ORDER,N,Z,SIGMA,SNRj,GAMMAv,KYO) 
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SNRj = SNR(2) 
NPRIME = gcr(SNRj,GAMMAV) 
SNRjbar = SNRj/NPRIME 
for r = 1:Order+2 
KY1(r) = gcumu(r,SIGMA,SNRjbar,GAMMAV) 
end 
fZvl = gdpz(ORDER,N,Z,SIGMA,SNRj,GAMMAv,KY1) 
plot(Z,fZinfO,'-',Z,fZvO,'--',Z,fZinfl1,'-',Z,fZv1,'--') 
grid 
xlabel(’Normalized Output Voltage, 2’) 
ylabel(’Probability Density, f(z)’) 


This program computes and plots the ROC data points for 
the Gaussian square law system. It is programmed to provide 
accuracy up to the 4th grouping of Edgeworth’s series. To 
obtain any other order of accuracy, the local variable 
ORDER which is set to 4 and the function file gdfz.m must 
be modified accordingly. 


oo o9 o? o? o? oe 


a? 


disp(’Choices:’) 
disp('1. Computing Pfa and Pd given a normalized threshold.’ ) 
disp(*2. Pd vs Pfa curves.’) 


disp(’ ”) 
CHOICE = input(’Select one of the above:  ') 
if CHOICE == 
ORDER = 4 
MORE = 1 
while MORE == 
N = input(’Number of samples integrated =  ') 
SNR = input(’Input SNR = ") 
SIGMA = input(’Input rms noise voltage =  ') 
GAMMA = input(’GAMMA = ') 
VT = input(’Normalized threshold voltage, VT = ^") 
disp(’Computing ... ’) 


for r = 1:ORDER+2 
KYO(r) = gcumu(r,SIGMA,0,GAMMA) 
end 
NPRIME = gcr(SNR,GAMMA) 
SNRjbar = SNRj/NPRIME 
for r = 1:ORDER+2 
KY1(r) = gcumu(r,SIGMA,SNRjbar ,GAMMA) 
end 
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Pfa = upfz(ORDER,N,VT,SIGMA,KYO,1) 


Pd = upfz(ORDER,N,VT,SIGMA,KY1,1) 
disp(’Do you wish to compute Pd vs Pfa for another VT?’) 
MORE = input('Type 1 if YES, 0 if NO:  ") 
end 
else 
N = input( ’Number of samples integrated = ') 
GAMMA = input( “GAMMA = 7) 
disp(’Enter up to a maximum of 3 SNR values as elements of 
Ei) 
SNR = input(’SNR [] = ^") 
SIGMA = input(’Input rms noise voltage = ') 


disp(’Enter range & resolution of z to be used in 
computing’) 
disp(’the ROC.’) 
range = input(’[zmin,zmax,zres]’) 
Z = range(1):range(3):range(2) 
for r = 1:ORDER+2 
KYO(r) = gcumu(r,SIGMA,0,GAMMA) 
end 
NPRIME = gcr(SNR, GAMMA ) 
SNRjbar = SNR ./NPRIME 
for r = 1:ORDER+2 
KY1(r) = gcumu(r,SIGMA,SNRjbar,GAMMA) 
end 
Pfa(j,:) = upfz(ORDER,N,Z,SIGMA,KYO,0) 
Pd(j,:) = upfz(ORDER,N,Z,SIGMA,KY1,0) 
if length(SNR) == 
semilogx(Pfa,Pd,’-’) 
elseif length(SNR) == 
semilogx(Pfa(1,:),Pd(1,:),Pfa(2,:),Pd(2,:)) 
elseif length(SNR) -- 3 


semilogx(Pfa(1,:),Pd(1,:),Pfa(2,:),Pd(2,:),Pfa(3,:),Pd(3,:)) 
end 
grid 
xlabel('Probability of False Alarm’) 
ylabel(’Probability of Detection’) 
end 


— — — — — — a cs cs i ws cs cs es es cs cs re ee es ee ee a 
—_ rm a PP rr mm ee s mm 
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de GCUMU.M FUNCTION 
function K - gcumu(r,SIGMA,SNR,GAMMA) 


% This function computes the rth cumulant of the output pdf 
% of a square law detector with Gaussian IF and video filters. 
% GAMMA is the ratio of the IF bandwidth to twice the video 
$ bandwidth. SIGMA is the input rms noise. 

% Reference: Thesis Eqn (10). 


G12 = 1 + 2*GAMMA“2 
(sqrt(G12) + 1)^r 
(sqrt(G12) - 1)°r 
BN1 = (2 * SIGMAS2) r 
rsl 


Q 
- 
NO 
a 
HM 


B = (BN1 * BN2)/BD1 

CN1 = BD1 

CD1 = G12PR + G12MR 

C = sqrt(CN1/CD1) 

K =B * (1 + r*SNR*G12”0.25) 


4. GDFZ.M FUNCTION 
function fz = gdfz(ORDER,N,Z,SIGMA,SNR,GAMMA,K) 


This function computes the probability density for specified 
values of the Gaussian square law system normalized output, 


K is an array comprising the cumulants of the pdf. 
Reference: Thesis Eqn (29). 


oe oe N 00 oe 


for j = 1:ORDER 
SUMn = 0.0 
SUMd = 0.0 
for m = 1:N 
SUMn = SUMn + K(m,j+2)/N 
SUMd = SUMd + K(m,j+2)/N 


N = SUMN/(N^(j/2) x SUMda 7727259 
en 

fz = zeros(1,length(Z) ) 

SUM = zeros(ORDER,length(Z)) 

X = (SIGMA*Z - N*K(1))/sqrt(N*K(2)) 
CONST - SIGMA/sqrt(N*K(2)) 
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PSIO = exp(-(Z .*2)/2)/sqrt(2*p1) 
for j = 1:ORDER 
if j = 1 
SUM(1,:) = CONST * PSIO .*( 1 + hmt(Z,3) * C(1)/6 ) 
elseif j == 
SUM(2,:) = CONST * PSIO .*( hmt(Z,4) * C(2)/24 + hmt(Z,6) 
AAA EN, 
elseif j == 
SUM(3,:) = CONST * PSIO .*( hmt(Z,5)*C(3)/12 + hmt(Z,7)* 
C(1)*C(2)/l144 + hmt(Z,9)*(C(1))^3/1296 ) 
elseif j == 4 
SUM(4,:) = CONST * PSIO .*( hmt(Z,6)*C(4)/720 + hmt(Z,8)* 
(ce?) 2/1152-4+ C( 1) *C@eay720 ) + hmt(Z,10)+* 


(CC T5) a? * C257 11 7 2:8 + 
hmt(Z,12)*((C(1))^4/31104) 
end 
fz = fz + SUM(j,:) 
end 
5% UPFZ.M FUNCTION 


function UPTE = upfz(ORDER,N,Z,SIGMA,ERROR) 


This function computes the Pfa or Pd of a Gaussian square 
law detector for a specified threshold normalized by the rms 


of oe 


noise voltage at the input to the detector. If the argument 
ERROR is set to 1, the function will return a two element 
row vector. The second element indicates the order of 
accuracy of the Pfa or Pd value given by the first element. 
If ERROR is set to O, then the function will return the 
the value of Pfa or Pd without any indication of the order 
of accuracy. 

Reference: Thesis Eqn (37) and Eqn (38). 


o9 o9 oo o? o? o9 o9 oe 


oo 


= 
for j = 1:ORDER 
SUMn = 0.0 
SUMd = 0.0 
for m = 1:N 
SUMn = SUMn + K(m,j+2)/N 
SUMd = SUMd + K(m,2)/N 
end 
k) = SUMn/(N^(j/2) * SUMA^((j+2)/2)) 
en 


ZZ = (SIGMA*Z - N*K(1))/sgrt(N*K(2)) 
PSIO = exp(-ZZ .*2)/sqrt(2*pi) 

2M = ZZ/sqrt(2) 

SUM1 = 0.5*erf(-INF,ZM, ‘high’ ) 
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TERMi = hmt(ZZ,2)*C(1)/6 

TERM2 = hmt(22,3)*C(2)/24 + hnt(22,5)*C(1)*2/72 

TERM3 = hmt(ZZ,4)*C(3)/120 + hmt(ZZ,6)*C(1)*C(2)/144 + 
hmt(ZZ,8)*C(1)*3/1296 

TERM4 = hmt(ZZ,5)*C(4)/120 + hmt(ZZ,7)*C(2)*2/1152 + 


hmt(ZZ,7)*C(1)*C(3)/720 + hmt(ZZ,9)*C(1)*2*C(2)/1728 
hmt(ZZ,11)*C(1)^4/31104 
SUM2 = PSIO .*(TERM1 + TERM2 + TERM3 + TERM4) 


if ERROR == 1 
UPTE(1) = 1 - (SUM1 - SUM2) 
UPTE(2) = PSIO .*TERM4 
else 
UPTE = 1 - (SUM1 - SUM2) 
end 
e GCR.M FUNCTION 


function GRHO - gcr(SNR,GAMMA) 


= á < es Ss Ss Ss SS nn nn ee in i PR 
KK KA A A KA KA sc ee ce RÁ A A eee ee 


oe 


This function computes the collapsing ratio caused by 
insufficient video bandwidth in a Gaussian square law 
system. GAMMA can be a vector but SNR must be single 
alued. 

Note that GRHO is also the equivalent number of independent 
Samples averaged by the video filter. 

Reference: Thesis Eqn (26). 


oe o9 oo «& ov o9 oo 


oo 


1.4 2*GAMMA .?2 

1 + GAMMA .%2 

B = 2*SNR*sqrt(B12 ./B11) - sqrt(B12) 
C = -2*SNR*sqrt(B12) 

GRHO = 0.5*(-b + sqrt(B .*%2 - 4*C) 


UJ 
FP 
- 
uod 


72 HMT.M FUNCTION 
function H = hmt(Z,m) 


This function computes the mth order Hermite polynomial for 
values of m less than or equal to 12. 


if m == 0 
H = 1 
elseif m == 
H = Z 
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elseif m == 2 
H-2.^2-1 
elseif m -- 3 
H = 2 .^3 - 3%Z 
elseif m == 4 
H = Z .*4 - 6%2 .^2 + 3 
elseif m == 5 
H= Z .°5 - 10%2 .^3 + 15%2 
elseif m == 6 
m 6 — 15)2 42 4 + 45*Z2 .*2 - 15 
elseif m == 7 
H = Z .^7 =- 21%Z .^5 + 105*Z .^3 = 105%2 
elseif m == 8 
IS AS SAA 6 1 210*57 .^4 - 420*7 .^2 + 105 
elseif m == 9 
H= Z .*9 - 36*%2 .*7 + 378%*Z .^5 - 1260*Z72 .^3 + 945*Z 
elseif m == 10 
H =Z .^10 - 457Z .^8 + 630ž*Z .^6 >= 3150*Z2 .^4 + 4725*Z .*2 
- 945 
elseif m == 11 
HEUS - 55*Z .^9 -FO990*7 .^7 — 6930*Z .^5 + 17325*Z .^3 
- 10395%2Z 
elseif m == 12 


MM A. 12 — GO*Z .^l10 «^ 1485*2 .^8 — 13860*2 .“6 
DAD *Z 5° 4 2462370*Z2 .^2umb»10395 
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